Theory of Plasmon-Enhanced Metal Photoluminescence 
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Metal photoluminescence (MPL) originates from radiative recombination of photoexcited core 
holes and conduction band electrons. In metal nanostructures, MPL is enhanced due to the sur- 
face plasmon local field effect. We identify another essential process in plasmon-assisted MPL — 
excitation of Auger plasmons by core holes — that hinders MPL from small nanostructures. We 
develop a microscopic theory of plasmon-enhanced MPL that incorporates both plasmon-assisted 
enhancement and suppression mechanisms and derive the enhancement factor for MPL quantum 
efficiency. Our numerical calculations of MPL from Au nanoparticles are in excellent agreement 
with the experiment. 



I. INTRODUCTION 



Since its discovery by Mooradian in 1969 [l|, photolu- 
minescence of noble metals and, later, metal nanostruc- 
tures has attracted increasing interest The un- 
derlying mechanism of metal photoluminescence (MPL) 
is radiative recombination of photoexcited d-band holes 
and upper-lying s-band electrons — a process strongly 
suppressed in the bulk by a plethora of nonradiative 
transitions in noble metals Q. In confined metal struc- 
tures, MPL is considerably more efficient due to plas- 
mon enhancement of radiative transitions Q. During 
the past decade, a highly efficient single-photon and 
multiphoton-absorption-induccd MPL was reported from 
various plasmonic systems inclu ding spherical nanoparti- 
cles mPl(iJ3-p i M HQ, [Ilia, nanorods llBlIS- 
03 M, Eflll IfT^nanoshells EH, nanowires 1171 1 181 . 
and bowtie antennas fl4j . This continuing interest in 
plasmon-enhanced MPL is fueled, in large part, by its 
promising applications, e.g., in imaging of blood vessels 
[Tol | and cancer cells [lg and optical recording [28| . 

In this Letter, we present a microscopic theory of MPL 
from plasmonic nanostructures. To highlight the issue at 
hand, consider radiative recombination of a d-holc and an 
s-elcctron inside a spherical NP of radius a much smaller 
than the radiation wavelength A. Since d-holes in noble 
metals are strongly localized, this process can be viewed 
as radiation of a dipole located at the d-hole position 
r. Radiative decay rate of such a dipole, T^, is enhanced 
relative to free-space radiative decay rate, 7^, by the local 
field factor Q 
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where e(w) and eo are, respectively, metal and outside 
medium dielectric functions. The factor ([T]) has a peak 
when radiation frequency uj is close to surface plasmon 
frequency; the latter is defined as the pole in NP dipole 
polarizability, ol\ = a 3 (e — eo)/(e + 2e ) [see Fig. QJa)]. 
Importantly, Eq. ([TJ depends neither on NP radius a nor 
on d-hole position r, thus implying universal MPL en- 
hancement within wide NP size range, a -C A. 

However, the experiment points to a more complicated 
picture. Systematic studies of MPL size dependence for 
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FIG. 1. (a) Plasmon enhancement of d-hole radiative recom- 
bination for Au NP in water (solid line) and normalized ab- 
sorption (dashed line), (b) Auger plasmon excitation rate by 
a d-hole localized at the NP center. 



spherical NPs 0, Q revealed a strong emission reduction, 
by a factor 3-5, for small NP with diameters D < 3 nm. 
Furthermore, MPL from larger NPs was reduced as well 
0, [H , and the brightest NPs observed were those with 
D ~ 60 nm While for large NPs, plasmon enhance- 
ment is expected to weaken due to the retardation ef- 
fects, its decrease for smaller NPs indicates that another 
mechanism, sensitive to NP size, is at work. Here, we 
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identify this mechanism as excitation of Auger plasmons 
by d-holes. 

Namely, a localized d-hole can undergo a highly effi- 
cient nonradiative recombination with an s-electron ac- 
companied by excitation of surface plasmons with high 
angular momenta /. This process is analogous to the tra- 
ditional d-hole Auger scattering to s-band with a Fermi 
sea electron-hole pair being excited except now collec- 
tive, rather than single-particle, excitations populate the 
final state. Importantly, the Auger plasmon excitation 
rate, T^, strongly depends on the proximity of d-hole po- 
sition r to the metal boundary. For example, excitation 
of high-/ Auger plasmons by a d-hole located near the 
center of a spherical NP (r/a <C 1) is suppressed by sym- 
metry and only the dipole (/ = 1) plasmon is excited 
with the rate |29l - l3ll | oc a~ 3 ; but for a d-hole near 
NP boundary (r/a ~ 1), high-/ plasmons are efficiently 
excited causing increase by several orders of magni- 
tude. In small NPs, where even dipole Auger plasmons 
are efficiently excited [29M3l| . the plasmon-assisted non- 
radiative decay rate F^ can exceed (on average) the bulk 
decay rate 7j r . This leads to the size- dependent reduc- 
tion of plasmon-cnhanced MPL. In fact, d-holc radiative 
recombination is quenched near NP boundary, and so 
the light emanates mainly from NP central region. Be- 
low we derive general expressions for MPL enhancement 
factor applicable to arbitrary metal nanostructures and 
present explicit formulas for a spherical NP. Our numer- 
ical calculations of MPL quantum efficiency (QE) for Au 
NP indicate that MPL from small NPs is significantly 
reduced and that brightest NPs are those with D ss 60 
nm, in excellent agreement with the experiment [H, Q . 

Before we proceed, let us compare classical and mi- 
croscopic descriptions of MPL enhancement mechanism. 
Within the classical approach [H, the radiating dipole 
comprised of a d-hole and an s-clcctron is enhanced by 
the surface plasmon local field. Within the microscopic 
approach [H , MPL enhancement is caused by the appear- 
ance of plasmon-mediated radiative decay channel: the 
d-hole recombines nonradiatively by exciting a surface 
plasmon [29l - l3i| which then decays radiatively, leaving 
a photon in the final state. Within Fermi golden rule, 
the photon emission rate is obtained by incorporating 
plasmons as intermediate states of the transition matrix 
element. Naturally, both approaches yield the same en- 
hancement factor ([1]) for MPL intensity. However, exci- 
tation of high-/ Auger plasmons with lifetime of ~10 fs 
[3^ by localized d-holcs provides a new nonradiative de- 
cay channel that can suppress MPL from small plasmonic 
structures. Below, we incorporate Auger plasmons exci- 
tation processes within a microscopic model for MPL. 



II. THEORY 

To specify the scope of our theory, let us elucidate 
the relevant processes within the standard three stage 
MPL model [2|, H[. Initially, the incident light excites 



a nonequilibrium d-hole population which, during the 
second stage, undergoes fast energy and momentum re- 
laxation due to electron-electron and electron-phonon 
scattering. During the third stage, those d-holes that re- 
laxed to states with momenta p < pp, where pf is the 
s-band Fermi momentum, recombine with upper-lying s- 
band electrons via vertical (in momentum space) tran- 
sitions [sec cartoon in Fig. QJa)]. The overall MPL QE 
is given by the product q = q e q r qd, where q e , q r and qd 
are, respectively, the bulk QEs for d-holc excitation, re- 
laxation, and radiative recombination stages. The latter 
two stages are dominated by fast nonradiative processes 
in noble metals @], including highly efficient (in the ab- 
sence of energy gap between d-band and s-band) d-hole 
Auger scattering, which quench d-hole momentum relax- 
ation and radiative recombination and result in a very 
low bulk MPL QE of [IH1 q ~ 1(T 10 . 

In plasmonic structures, QEs for individual MPL 
stages are independently modified. For excitation stage, 
a change in q e can be substantial if some spectral overlap 
between the initial excitation and surface plasmon bands 
is present. Here we assume, however, that the incident 
light frequency is either well abov e ( for single-photon- 
absorption- induced MPLH-SHtHIlHil) or wel1 below 

(fo r multiphoton-absorption-induced MPL 0, 191-H21. lT3 - 
EH H3) the plasmon band so that q e remains unaffected. 
The QE for fast relaxation stage, q r , is not directly af- 
fected by plasmons although the increased d-hole surface 
scattering may lead to faster momentum relaxation in 
small nanostructures [33j]. However, q r is significantly 
enhanced for multiphoton (vs. single-photon) excitation. 
Indeed, the multiphoton d-hole excitation in noble metals 
is dominated b y th e following incoherent sequential pro- 
cess [ill, [2^, [22|.|27( : an incident photon creates a vacancy 
in the s-band below the Fermi level that is filled by a d- 
band electron promoted by a subsequent photon. In this 
way, d-holcs arc created with momenta predominantly 
below pp and hence need not undergo large momentum 
transfer before radiative recombination. In the absence 
of phase space restriction, q r drastically increases result- 
ing in a highly efficient multiphoton-absorption-induced 
MPL @, fflU, EHH El- A quantitative description of 
the above MPL dynamics is beyond the scope of our pa- 
per. Here, we focus solely on the third (recombination) 
MPL stage, characterized by quantum efficiency Qd that 
is strongly influenced by plasmons, and derive the QE 
enhancement factor, M = Qd/ qd- 

We start with excitation of Auger plasmons by a local- 
ized d-hole in a spherical metal NP of radius a. The decay 
rate of a photoexcitcd d-hole is given by Yd = 2Im£<j, 
where is d-holc self-energy due to electron-electron 
and electron-phonon interactions as well as interaction 
with light. In plasmonic nanostructures, the interac- 
tions between electrons (and holes) are described by 
dynamically-screened Coulomb potential U(u; r, r') = 
U (r -r') + 6U(uj;t,t'), where U = e 2 / (e(w)|r - r'|) 
is the bulk Coulomb potential (e is the electron charge) 
and SU contains contribution from surface excitations. 
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In spherical NP, SU takes the form [291-43 ll 
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where <pi m (r) = r Yj m (f ) (Yj m are spherical harmonics) 
and the factor 1 /e; = (2i + l)/[Ze+(i+l)eo] — l/e contains 
a pole corresponding to surface plasmon with angular 
momentum i. The full decay rate Td can be split into 
sum Yd = Yd T + + T^, where 7^ r is bulk nonradiative 
decay rate, T d is plasmon-cnhanccd radiative decay rate, 
and T d is plasmon-assisted nonradiative decay rate. The 
latter is determined by the contribution to d-hole self- 
energy due to intcrband transitions under potential SU, 
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where Mj^, = e(s|<^ m (r)|d) is the interband transition 
matrix element. To evaluate Mf£, we adopt the tight- 
binding approach used for calculation of local field cor- 
rections to the dielectric function |33436j. First, we ex- 
pand the d-holc Bloch function in the Wannier basis as 



'(Pdir-Tj), where tp d (r-Tj) 



the localized hole wave-function for atomic position rj 
and V is unit cell volume. Note that, for a given d- 
hole position, only one term effectively contributes to the 
sum. In this basis, the matrix element of a smooth (on 
atomic scale) function </>j m (r) between d-band and s-band 
states is e(s\<j>i m (r)\dj) = d sd ■ V<& TO (rj), where d sd = 
e(s | (r — Tj)\dj) is the intcrband dipole matrix element 
|3#j36j . The next step is to average (\d s d ■ V</>/ m (r)| 2 ) 
over orientations of d s d assuming isotropic angular dis- 
tribution: (d^ d d" sd ) = \d 2 sd 8^ y . Finally, using the relation 
27TT £ m |V^ m (r)| 2 = 1(21 + l)r 21 - 2 , we obtain 
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The surface-induced decay rate T d = 2ImE^ peaks at 
frequencies satisfying e(ui) + (1 + l/i)eo = indicating 
Auger plasmon excitation. A d-hole at the NP center 
(r = 0) decays into dipole plasmon only with the rate [29l — 
= -2 (d 2 d /a 3 ) Im[l/ei(w)]. In Fig.QJb), we plot 
to illustrate the huge (up to 10 5 ) difference between 
plasmon-assisted radiative and nonradiative decay rates. 

Turning to MPL QE, the radiated energy density S(uA 
of a dipole with excitation energy E sd is given by |38l - t40l | 



S(lo) oc 
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where Ed is d-hole full self-energy and Qd = F d /Td is QE 
for d-hole radiative recombination. Summing up Eq. ([5]) 
over empty d-band and occupied s-band states and re- 
turning to the momentum space j3^ - l36| . we obtain 



S(cj) (X uQ d (u)x'd(u) 



(6) 



where 



Xd(w) = —2 
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sd 



is intcrband susceptibility modified by the surface, 
E* — E d is vertical transition energy between d-band 
and s-band with dispersions E d and E*, respectively, 
/x = imEp d d s d/e is the interband momentum matrix ele- 
ment (m is electron mass), and / is the Fermi function. 
Importantly, the self-energy E^ in Eq. (O depends on 
w rather than on E^f [2iM3ll ]. which allows us to factor 
out Qd(w) in Eq. (|6]). The bulk radiated energy den- 
sity is recovered from Eqs. (|6j and ([7]) by replacing E^ 
with i-id = i(jnr + 7r)> yielding S cx uq d x'd(uj), where 
Id = Id/ {id + lY) an( i Xd(w) are, respectively, bulk 
QE and interband susceptibility. Note that in Eq. (0, 
the plasmon contribution to Ed plays no significant role 
due to the integration over off-resonant p-states, so that 
Xd(w) in Eq. ((BJ) can be replaced by its bulk counterpart 
Xd(w). The local enhancement factor, A(lj,t) = S/So, 
then takes the form 
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K ' q d q d {P + R-l) 



(8) 



where R(u),r) — T r d /j d and P(w,r) = T d /^ d both, in 
general, depend on system geometry and d-hole position. 

The average enhancement factor, M, is obtained by av- 
eraging the local enhancement factor A(u),r) over d-hole 
positions inside a nanostructure. Here we note that d- 
holes spatial distribution stays largely unaffected during 
the fast relaxation stage, so we can adopt the normalized 
distribution function F(wj,r) = 5(w,,r)/ J drS, where 
S(u)i,r) is the absorbed power density, 



S(u H ,r) = ^-e"(uj i )\E(u; i ,r)\ 2 , 



(9) 



and E(cj,,r) is the local field at incident light frequency 
LUi. Then the enhancement factor per radiating d-hole 
is given by the weighted average of A(ui,r) over metal 
volume 



M(uj, 



dVF{ui,T)A{w, 



(10) 



Equations (j5j)- (fTU)) represent our model of plasmon- 
assisted MPL that applies to any metal structure. Below 
we consider MPL from a spherical NP integrated over 
emission angles. Since the excitation and radiation stages 
are separated by the relaxation stage, the incident and 
emitted light polarizations are uncorrelated. In this case, 
R and P in Eq. (J5J depend only on radial coordinate r, 
and F(uJi,r) in Eq. (fTOj) can be replaced by its average, 
F(ui,r) = J dQF(uJi,r). Using Mie theory formulae 
for the electric fields inside a metal sphere [4JJ, we have 
F(ui,r) = F (u}i)R(u}i,r) with F Q = e"ki/e Q C ab , where 
C a b is NP absorption crossscction, ki = ujiy/eo/c is the 
wavevector and c is speed of light, yielding 
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FIG. 2. Surface-induced nonradiative decay rate (a) and local 
enhancement factor (b) near plasmon resonance (A = 550 nm) 
is shown vs. d-hole radial position in Au NP. 
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FIG. 3. MPL enhancement factor for Au NP with diameters 
in the range 2-40 nm (a) and 60-160 nm (b). The maximal 
enhancement is achieved for D — 60 nm NP. 



e"(^)h r dVR(ui,r)R(u>,r) 
C ab (u>i)e J q d [P(uj,r)+R(u>,r)-l} + l' [ ' 

For small NP (o < A), fl and P are derived, respectively, 
from Eqs. (JTJ and Q using Yd = 4d^w 3 v /eo/3c 3 , and 
Cab = 4:Trkia'{; the formulas for NPs of arbitrary size 
[4ll |4^ | are provided in the Supporting Information. 



III. NUMERICAL RESULTS 

Calculations were performed for an Au NP in water 
using standard Au parameters [37j along with bulk d-hole 
decay time of 40 fs @, d, rendering qd in the interval 
(3-4) xlO -7 within the plasmon band range 500-600 nm. 
Angular momenta up to l m ax = 100 were included in 
the evaluation of 1^ (in small NP, l max is restricted by 
app), and experimental Au dielectric function was used 
throughout. 

In Fig. [U[a), we plot normalized to 7^ against d- 
hole position r for ui near the dipole plasmon resonance 
(compare to Fig. [T]). For d-holc at the NP center, only 



the dipole Auger plasmon is excited, but for r ~ a, T p d 
increases by several orders of magnitude due to the high-Z 
contribution [sec Eq. (dJ)]. Near NP boundary, T p d exceeds 
significantly the bulk d-holc decay rate Y d lr , leading to 
drastic reduction of local enhancement factor A(lu, r) [see 
Fig. [UJb)]. This effect is especially strong in small NP 
(D = 2 nm) where and j d r are largely comparable 
in most of NP volume (i.e., away from NP center and 
its boundary). Note that for D smaller than the skin 
depth (« 20 nm for Au), d- holes are distributed nearly 
uniformly throughout NP volume and their distribution 
function is simply F = 3/47ta 3 . In this case, it is NP 
central region that mainly contributes to MPL, i.e., NP 
has a "bright" core and "dark" edge that expands toward 
NP center as its size decreases. 

In Fig. [3] we plot the average enhancement factor M, 
calculated from Eq. (fTT|) . for Au NPs with diameters 
between 2 nm and 160 nm. In a wide size range of 
4-100 nm, M varies only weakly (within factor of 2), 
and the strongest enhancement is achieved for D = 60 
nm NP, consistent with experiment For small NP 
(D = 2 nm), M falls by a factor ^4 relative to maxi- 
mum enhancement, also consistent with the data [1, Q. 



5 



Note that this remarkable agreement with experiment is 
achieved without use of fitting parameters. For larger 
NPs (D = 80 and 100 nm), the MPL peak is rcdshifted 
0, and with a further size increase (D = 160 nm), the 
plasmon band nearly disappears and MPL is diminished. 
Thus, the emergence of optimal NP sizes in the range of 
40-80 nm is due to the interplay between two different 
damping mechanisms: weakening of plasmon resonance 
in larger NPs due to electromagnetic effects (retardation) 
and MPL quenching in smaller NPs due to quantum ef- 
fects (excitation of Auger plasmons) . 



with 

Df = eji(xi) [x 2 hi(x 2 )}' - e hi(x 2 ) [xiji(xi)]' , 

D? = ji(xi) [x 2 hi(x 2 )}' - hi(x 2 ) [xijjfxi)]' > ( 15 ) 

where hi (x) = ji (x) + iyi (x) is the spherical Hankcl func- 
tion and we adopted the notations x\ = (w^e/c) a and 
x 2 = (w^/eo/c) a (prime stands for derivative). Nor- 
malizing Td by radiative decay rate in outside medium, 
Yd = id 2 uj 3 y/eo/3c 3 , and averaging over dipole orienta- 
tions, we obtain 



IV. CONCLUSION 

In summary, we presented a microscopic theory for 
plasmon-enhanced metal photoluminescence. We have 
shown that, in small nanostructures, plasmon-enhanced 
luminescence is hindered due to highly efficient excita- 
tion of Auger plasmons by photoexcited d-holes. Our 
numerical calculations of the MPL enhancement factor 
for spherical Au nanoparticles are in excellent agreement 
with experiment. 

This work was supported in part by NSF Grants DMR- 
1206975 and HRD-0833178. 



SUPPORTING INFORMATION 

Here we present expressions for the surface-induced 
decay of an excited dipole comprised of a localized d- 
hole and s-electron inside a spherical NP of radius a. 
This problem is mathematically similar to the decay of a 
dipole placed inside a dielectric sphere [IH, so we adopt 
that solution here while noting the differences between 
dielectric and metal spheres. The decay rate is related to 
the electric field Green dyadic, G(w;r,r'), at the dipole 
position r as j38j 

Pi = ^Imd.G(aj;r,r) • d, (12) 

cr 

where d = dw is the dipole moment and averaging over its 
orientations n is implied (the effects of electron-electron 
and electron-phonon scattering can be incorporated by 
adding the bulk decay rate j d r , sec main text). The 
Green dyadic has free and scattered parts, G = G° + G s , 



where G°„(cj; i^r) 
has the form [42( 



5^k/6Tr (here k = uiy/e/c), and G s 
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EiKl + l) 



3i{x) 



rr 



Ei 



+ W(x) 



8G + <p<t> 



, (13) 



where (r, 0, cf>) are unit vectors in spherical coordinates, 
ji{x) is the spherical Bessel function and x = kr. The 
Mie coefficients E[ and Ft are given by 

Ei = [e hi(x 2 ) [xihi{xx)}' - ehi(x\) [x 2 hi(x 2 )}'] /Df, 

F = [hi(x 2 ) [xxhixx)] 1 - hi{xx) [x 2 hi(x 2 )}'] /Df 1 , (14) 
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Id V £ o 

+ E, 
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(16) 



The radiative decay rate T r d is derived by integrating the 
Poynting vector over a sphere with infinite radius [42[ , 
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(17) 



with coefficients Df and Df 1 given by Eq. (fT5")) . Note 
that for dielectric sphere (e > 0), the two rates (fTB]) and 
(JT7J) coincide [13] . For metal sphere with complex dielec- 
tric function e(uj) considered here, radiative rate is still 
given by Eq. (|T7| , while now = + T^. In our nu- 
merical calculations, we used the above normalized rates, 
R(ui,r) = T d /Y and P(w,r) = T d /Y- I n the longwave 
limit, ak <C 1, it can be readily checked that those coin- 
cide with the corresponding expressions in the main text. 

Finally, the NP absorption crossection can be found 
as the difference between the extinction and scattering 



crossections, C a b 
way [il| 

2tt 



C sc , defined in the standard 



C ex = -jReJ2(U + l)(ai+bi) 



2tt, 



i 



C sc ^^RcY,(2l + l)(\ai\ 2 + \bi\ 2 
fc o i 



where Ai and Bi are the Mie coefficients, 



(18) 



A i = [ e Ji(xi) [x 2 ji(x 2 )}' - e ji(x 2 ) [xiji(xi)Y] /Df, 
Bi = [ji(xi) [x 2 ji(x 2 )Y - ji(x 2 ) [xdMY] /DY 1 . (19) 

These expressions were used for the calculation of d-holc 
distribution function. 
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